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A partial wave analysis of XcO — > in ip{2S) — > 7XcO decay is presented using a sample 

of 14 million tp{2S) events accumulated by the BES II detector. The data are fitted to the sum of rel- 
ativistic covariant tensor amplitudes for intermediate resonant decay modes. From the fit, significant 
contributions to XcO decays from the channels /■o(980)/o(980), /o(980)/o(2200), /o(1370)/o(1710), 
K* (892)° K* (892)°, J^q* (1430)^-0* (1430), {U30)K^ (1430) + c.c, and Ki{1270)K are found. 
Flavor-SU(3)-violating 7^1(1270) — 1400) asymmetry is observed. Values obtained for the masses 
and widths of the resonances /o(1710), /o(2200), /o(1370), and -K'o(1430) are presented. 



PACS numbers: 13.25.Gv, 12.38.Qk, 14.40.Gx 
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I. INTRODUCTION 

Exclusive heavy quarkonium decays constitute an important laboratory for investigating perturbative QCD. Com- 
pared to J/V' decays, relatively little is known concerning Xcj (J ~ 0,1,2) decays More experimental data on 
exclusive decays of P-wave charmonia are important for a better understanding of the nature of Xcj states, as well 
as testing QCD based calculations. Further, the decays of Xcj, in particular and Xc2, provide a direct window on 
glueball dynamics in the 0^"*" and 2++ channels, as the Xcj hadronic decays may proceed via cc — > — *■ gggg 01 ■ 

Amplitude analysis of XcJ decays is an excellent tool for studying charmonium decay dynamics. Knowledge of the 
quantum mechanical decay amplitude allows one to investigate not only the intermediate resonant decay modes but 
also to properly account for the interference effects between different resonances. 

In this paper, partial wave analysis results of Xco — ^ 'n^'K~ K'^K~ in ipi'^S) — > 7XcO decays using 14 million ^p{2S) 
events accumulated at the BESII detector are presented. In previous studies of this channel, only the decay modes 
Xco K+K*{892)°n- + c.c. and K*{892fK*{892)° were measured Il|3|. Here additional information from partial 
wave analysis is very important. 

II. BBS DETECTOR 

BES II is a large solid-angle magnetic spectrometer that is described in detail in Ref . . Charged particle momenta 
are determined with a resolution of cTp/p ~ 1.78%-\/l + (p in GeV/c) in a 40-layer cylindrical drift chamber. Particle 
identification (PID) is accomplished by specific ionization (dE/dx) measurements in the drift chamber and time-of- 
flight (TOF) measurements in a barrel-like array of 48 scintillation counters. The dE/dx resolution is (JciE/dx = 8.0%; 
the TOF resolution is cttof — 180 ps for Bhabha events. Outside of the TOF system is a 12 radiation length lead-gas 
barrel shower counter (BSC), operating in self-quenching streamer mode, that measures the positions and energies 
of electrons and photons over 80% of the total sohd angle. The energy resolution is as/E = 22%/\/E {E in GeV). 
Surrounding the BSC is a solenoid magnet that provides a 0.4 T magnetic field in the central tracking region of the 
detector. Three double-layer muon counters instrument the magnet flux return and serve to identify muons with 
momentum greater than 0.5 GeV/c. They cover 68% of the total solid angle. 

In this analysis, a GEANT3 based Monte Carlo (MC) simulation package (SIMBES) with detailed consideration 
of detector performance (such as dead electronic channels) is used. The consistency between data and MC has been 
checked in many high purity physics channels, and the agreement is quite reasonable 

III. EVENT SELECTION 

The selection criteria described below are similar to those used in previous BES analyses HQ. 

A. Photon identification 

A neutral cluster is considered to be a photon candidate when the angle between the nearest charged track and 
the cluster is greater than 15°, the first hit is in the beginning six radiation lengths, and the difference between the 
angle of the cluster development direction in the BSC and the photon emission direction is less than 30°. The photon 
candidate with the largest energy deposit in the BSC is treated as the photon radiated from the ip{2S) and used in a 
four-constraint (4-C) kinematic fit to the hypothesis ip{2S) 'yir^ ■k~ K~ . 

B. Charged particle identification 

Each charged track, reconstructed using MDC information, is required to be well fit to a three-dimensional helix, 
be in the polar angle region | cos^mdcI < 0.80, and have the point of closest approach of the track to the beam axis 
be within 2 cm of the beam axis and within 20 cm from the center of the interaction region along the beam line. 
For each track, the TOF and dE/dx measurements are used to calculate x^ values and the corresponding confidence 
levels for the hypotheses that the particle is a pion, kaon, or proton. 
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C. 



Event selection criteria 



Candidate events are required to satisfy the following selection criteria: 

(1) The number of charged tracks is required to be four with net charge zero. 

(2) The sum of the momenta of the two lowest momentum tracks with opposite charges is required to be greater 
than 650 McV/c; this removes contamination from ■0(25') tt^tt^ J/ip events. 

(3) The confidence level for the 4-C kinematic fit to the decay hypothesis ip{2S) ^tt^tt^ is required to be 
greater than 0.01. 

The combined confidence level determined from the 4-C kinematic fit and PID information is used to separate 
77r+7r~7r+7r~, ^K'^ K~ K'^ K~ , and the different possible particle assignments for the ^tt^tt^ final states. This 
combined confidence level is defined as 



where f{z;ndfau) is the probability density function, Xa/i the sum of the values from the 4-C kinematic fit 
and those of the four track PID assignments, and ndfaii is the corresponding total number of degrees of freedom. 
For an event to be selected, the combined confidence level of -fTr~^TT~ K^K~ must be larger than those of the other 
possibilities. In addition, the PID confidence level of each charged track must be > 0.01. 

Further rejection against Kg — > tt+tt" is obtained by requiring that any Tr~^TT~ K^K~ combination with M^^+tj- in 
the interval (497 ± 50) MeV/c^ should have r^y < 5 mm, where r^y is the distance from the beam axis to the tt+tt" 
vertex. 

The invariant mass distribution for the tt~^'!t~ K~ events that survive all the above selection requirements is 
shown in Fig. 1. There are clear peaks corresponding to the Xcj states. The highest mass peak corresponds to ip{2S) 
decays to four charged track final states that are kinematically fitted with an unassociated, low energy photon. 



FIG. 1: The tt'^-k K^K invariant mass spectrum. There are three clear Xcj peaks. The highest mass peak corresponds to 
'0(25') decays to four charged tracks final states that are kinematically fitted with an unassociated, low energy photon. 

An additional five-constraint (5-C) kinematic fit is made with the invariant mass of the tt^t:^ being con- 
strained to the Xco mass. After requiring the confidence level for the 5-C fit to be greater than 0.01, the tt^tt~ K'^K~' 
invariant mass distributions for data and MC, shown in Fig. 2, are obtained, where the tt^t:~ K'^ K~ invariant mass 
obtained from the 4-C fit and the generated mass and width of the XcO are fixed to PDG values 0. The agreement 
found in this comparison indicates a clean data sample suitable for partial wave analysis in which the four-momentum 
information from the 5-C fit will be used. Finally, 1371 ip{2S) JXcO, XcO ^ tt~^tt~ K^K~ candidate events are 
selected after all the above criteria. 

Figure 3 shows the individual K^K^, tt+tt^, Ktt, Ktttt, and KKtt invariant mass distributions. There are two 
strong symmetric structures at about 1.75 GeV/c^ and 2.2 GeV/c^ and some evidence for /o(980) in the K^K~ mass 
distribution, and there are clear p{77Q) and /o(980) peaks and a smaller one at about 1.3 GeV/c^ in the TT~^Tr~ mass 
distribution. There is no obvious structure for the KKt: mass spectrum in Fig. 3. 

Fig. 4(a) shows the scatter plot of K^K~ versus tt+tt" invariant mass for selected ip{2S) 'yXcO, XcO — > 
'!T~^Tr~ K~ events, which provides further information on the intermediate resonant decay modes for 
{tt^tt^){K^K^) decay. For instance, it can be seen that the /o(980) tt+tt^ mainly couples with /o(980), /o(1710), 
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FIG. 2: Tv^TT K'^K invariant mass distributions from the 4-C fit after requiring Probsc > 0.01, where the error points are 
data and the histogram is MC. 
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FIG. 3: The individual (a) K'^K , (b) tt^tt , (c) Ktt, (d) Ktvtt, and (e) KKtt invariant mass distributions after application 
of all selection criteria. 



and /o(2200), which decay to K^K~, and the concentration of events at large K^K^ mass is also associated with 
/9(770) decays to tt+tt^. 

Figure 4(b) shows the scatter plot of A'+7r~ versus A'~7r+ invariant masses. There are obvious clusters due to 
Jsr*(892)° and around 1.43 GeV/c^, where there are several known resonances, and evidence for [Kn) structures at 
about 1.7 GeV/c^ and 1.95 GeV/c^, as well as an enhancement above 2 GeV/c^. There is also a possible small 
accumulation outside the /\*(892)°A:*(892)" cluster, which may be due to the broad S-wave structure k. 

The K-KTT mass distributions for events where one i^T^Tr^ combination is in the K*{%^2) mass region 896±60 MeV/c^ 
and where the tt+tt" mass in the /c(770) mass range from 700 to 850 MeV/c^ are shown in Fig. 5. Strong i4ri(1270) 
signals are observed in both cases, and there is also a weak peak around 1.4 GeV/c^ for the A'*7r decay mode. 



IV. ANALYSIS METHOD 



We have carried out a partial wave analysis using relativistic covariant tensor amplitudes constructed from Lorentz- 
invariant combinations of the four-vectors and the photon polarization for ipi'^S) initial states with helicity ±1 
For ip{2S) 7XcO, XcO ^ tt^tt~ K~^K~ , the general form for the decay amplitude is 
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FIG. 4: The scatter plots of (a) K^K versus tt+tt and (b) K'^n versus K 7r+ invariant mass for selected tp{2S) — > 
IXcO, XcO — »■ 7r+7r~if+A'~ events. 
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FIG. 5: K-rvn mass distributions for events where (a) one K^n^ combination is in the mass region 896 ± 60 MeV/c^ and (b) 
the TT^TT^ mass is in the range 700 - 850 MeV/c^. 



where ^^{mi) is the ip{2S) polarization four-vector, 6^(7712) is the polarization four-vector of the photon, and U^'^ is 
the partial wave amplitude with coupling strength determined by a complex parameter A^. For the photon polarization 
four-vector e^, with photon momentum q, there is the usual Lorentz orthogonality condition e^q'^ = 0. We assume the 
Coulomb gauge in the 4'{2S) rest system with momentum p^/,, i.e., e^p.^ ~ 0. Then we have 

E*/ N / X , quKyi+K^qi,, K K 
e^(m2)e^-(m2) = -g^,^' + — (q- K)^ '^'''^''' = 

7712 \'d J 

with K = p.^ — q and e^K'^ = 0. We know that 

2 

mi —1 

so the radiative decay cross section is: 

2 2 

d$ " 2 



mi — 1 m2 — 1 



2 ■ 



where 



p. . _ p* — A - A* 
1 2 

Ai=r 
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The partial wave amplitudes Ui for the intermediate states used in the analysis, such as the K* (892)'^ K* (892)° are 
constructed with the four-momenta of the tt^ , 7r~ , , and K~ , and their specific expressions are given in Ref. 0| . 
For an intermediate resonance, the corresponding Breit-Wigner propagator is denoted by a function: 

1 



BW = 



AP 



iMT ' 



where s is the invariant mass-squared and Af , F are the resonance mass and width. Angular momenta L up to 2 in 
the production process are needed, but higher L give negligible contributions. Standard Blatt-Wcisskopf centrifugal 
barrier factors [HSl Etre included using a radius of interaction of 0.8 fm, though results are insensitive to this radius. 

The relative magnitudes and phases of the amplitudes are determined by an unbinned maximum likelihood fit. 
The basis of likelihood fitting is calculating the probability that a hypothesized probability distribution function 
(PDF) would produce the data set under consideration. If the probability to produce event i, characterized by the 
measurements Xi, is P{xi), then the joint probability density for observing the N events in the data sample is 

N 



C = Y[Pia 



1=1 



The normalization condition for P{xi) is that its integral over its domain must not depend on the values of the fit 
parameters. Suppose the differential cross section ( — ) is the unnormalized PDF for producing event i. Then 



da 



where the integration is over the domain of 



For the purpose of normalizing the PDF used in fitting the experimental data set, several million MC 'ip{2S) 
7XcO, XcO TT+TT^/i +/i ^ phase space events were generated. The events undergo detector simulation and are passed 
through the same analysis procedure as are the experimental data, and thus the distribution of the MC events passing 
into the final stage of analysis contains the acceptance information. The normalization integral is then computed as: 



1 



Nmc 



Nmc / , 

El da 
\ 7^ 




where Nmc is the number of accepted MC events. 

Background events obtained from MC simulation are included in the fit, but with the opposite sign of log likelihood 
compared to data. These events are used to cancel the backgrounds within the data sample in the maximum likelihood 
fit. This technique of background treatment has been used in analyses of Crystal Barrel data (see Ref. ^^^^ several 
previous BESII publications [Hill [11. Then finally we have the logarithm of C 



N 



da 



N 



El- 

i=l 



^ Pjk ■ Fjk 



Nmc 



N, 



MC V w 



i'=l ^ j,k 



where N is the total number of data and background events. 

For the '4^{2S) — > 7XcO, XcO ^ ti^t:^ K^K~ process, small backgrounds remaining arise mainly from ip{2S) 



"TT , ■0(25') 7XcO, XcO 



'TT TT+TT , i!{2S) JXcl, Xcl 



K+K- 



-K+K-, i;{2S) 

and '0(25') 7XcOi XcO ^ KgK^Ti^ . The number of background events is estimated to be about 29, only a 
few percent of the data sample, from detailed exclusive and inclusive MC simulations, and therefore the effect of 
background is expected to be minor in the partial wave analysis. In the MC simulation we increase the amount of 
background by a factor of 10 and then multiply by a normalization factor of 0.1 to reduce the statistical fluctuation 
of the background events. 

The optimization of the free parameters within the amplitude is done using FUMILI , which also gives the 
fitting error matrix. Technically, rather than maximizing ln£, this package minimizes 5 = — ln£. In the minimizing 
procedure, a change in log likelihood of 0.5 represents a one standard deviation effect. 
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V. ANALYSIS RESULTS 



Table I shows the decay modes considered in the partial wave analysis, which are motivated by the structures seen 
in the scatter plots of Figs. 4(a) and (b) and in projections of Figs. 3 and 5, changes AS" in log likelihood when 
the component is dropped from the fit, and the corresponding statistical significances. The partial wave amplitude 
improves log likelihood by more than 5 in most cases. The decay modes /o(980)/o(1710), /2(1270)/2(1270), and 
/o(1710)/o(1370), where for the latter /o(1710) decays to tt+tt" and /o(1370) decays to K^K^ , are shown for 
completeness, but they are not very significant. The significances are calculated from comparing the difference 
between the S values of the fits with and without the component. 



TABLE I: Decay modes fitted in the partial wave analysis, changes AS in log likelihood when the component is dropped from 
the fit, and the corresponding statistical significances. The errors are statistical only. 



Decay mode 


Fitted events 




Significance 


{■K TT )(K ]\ ) 










z/ .y zt 0. ( 


1 p; 7 

10. ( 


O.OCT 




1/1 7-1-7 


o.z 


Z.ocr 




77.1 ± 13.0 


27.3 


7.1cr 


/o(1370)/o(980) 


26.9 ± 10.0 


14.6 


5.O0- 


/o(1370)/o(17f0) 


60.6 ± 15.7 


23.5 


6.5(7 


/2(1270)/2(1270) 


5.9 ±4.1 


5.8 


3.0a 


a/o(1710) 


46.7 ± 13.4 


22.2 


6.3(7 


f7/o(2200) 


23.9 ± 8.8 


8.5 


3.7a 


/o(1710)/o(1370) 


4.6 ±4.9 


2.5 


1.7a 










if* (892)" A-* (892)" 


64.5 ± 13.5 


31.1 


7.1a 


7^* (892)° A'* (1680)" + C.C. 


40.5 ± 13.3 


21.0 


5.6a 


KS{U30)KS (14:30) 


82.9 ± 18.8 


28.0 


7.2a 


A'2*(1430)J^2*(1430) 


9.2 ±5.3 


7.1 


3.3a 


A-Q* (1430)J^2* (1430) +C.C. 


62.0 ± 12.1 


40.6 


8.7a 


A-Q* (1430)^-0* (1950) +C.C. 


71.0 ± 19.1 


22.7 


6.4a 




106.8 ± 16.7 


39.2 


8.6a 


K* {892f K* {2S00f + c.c. 


115.7 ± 19.4 


45.2 


8.8a 










Ki{1270)'^K- +C.C. 


153.0 ± 19.5 


102.2 


13.2a 


Ki{UOO)+ K- + c.c. 


19.7 ±8.9 


6.9 


2.7a 


K{U60)+K- +C.C. 


79.7 ± 16.8 


39.3 


8.2a 



The statistical uncertainties on the measurements shown in Table I arc derived from the uncertainties in the fitting 
parameters. Recall that the parameters used in the fitting function arc magnitudes and phases of the various processes. 
The numbers of events arc derived from the parameters using numerical integration of the individual amplitudes. Using 
the numerical expression of the total cross section a defined in section IV, the number of fitted events, Ni for an 
intermediate decay which has one partial wave amplitude Ui, is given by the integral over phase space (i.e. the sum 
over the MC events) of the PDF for that decay: 

H Nmc / \ 

= '-^ N , 

a 

where N is the number of events after background subtraction. For the intermediate decay with two partial wave 
amplitudes Ui, Uj, we replace Pa ■ Fa with Pa ■ Fa + Pij ■ Fij + Pji ■ Fji + Pjj ■ Fjj, where Pij ■ Fij and Pji ■ Fji 
correspond to the interference terms; and those with more amplitudes can be deduced similarly. 

A numerical scheme is used to extract the uncertainties on the numbers of fitted events. The standard deviations 
and the covariances of the fit parameters are obtained from FUMILI. Next, taking into account the correlation 
coefficients, a thousand sets of fit parameters are generated using a random number generator CORGEN These 
random numbers arc distributed as correlated Gaussian distributions. The number of fitted events for each process is 
then calculated with each of the thousand sets of parameters, and histogrammed. A Gaussian function is then fit to 
each of these histograms, and the standard deviations are determined. 



8 



The mass projections in K^K 
between the data and the fit. 



Kt:, and Kt^t: are shown in Fig. 6. There is a reasonable agreement 
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FIG. 6: Mass projections on (a) K+ K~ , (b) tt+tt^, (c) Kti and (d) Ktitt for 7/'(25') 
represent the fit result, and the points with error bars are data. 



XcO, XcO 



■K^-R K . The histograms 



A. XcO ^ {'K+-K-){K+K-) 

We begin by discussing the XcO (7r"*'7r~)(i<'+if ~) decay modes. In Fig. 6(b), the p(770) due to the decay of 
A'i(1270) AX770) and a strong /o(980) are observed. Along the /o(980) (-^ tt+tt-) band in Fig. 4(a), there are 
several enhancements, one at the K~^K~ threshold and others around 1.75 GeV/c^ and 2.2 GeV/c^, which correspond 
to the /o(980)/o(980), /o(980)/o(1710), and /o(980)/o(2200) modes listed in Table I. 

In the fit, the /o(980) is described with the usual Flatte formula the parameters used are those of 

Ref. ^|. At higher tt+tt" mass, a weak signal around 1.3 GeV/c^ is visible. It is mainly from couplings of the 
/o(1370) with /o(980) and /o(1710) which decay to K+R-. A free fit to /o(1370) gives a fitted mass of 1265 ± 30 
MeV/c^ and a width of 350 ± 100 McV/c^ with large statistical errors. Apart from the above structures, a correlation 
between the low tt+tt" mass enhancement and K~^K~ mass greater than 1.5 GeV/c^ also appears in Fig. 4(a). We 
describe it by adding cr/o(1710) and cr/o(2200) decay modes, where the parameterization of a J/ip —f wtt+tt" analysis 
(Adler zero parameterization) ^2] is adopted to describe the broad S-wave a. The cr/o(1710) and /o(1370)/o(1710) 
decays dominate the production of /o(1710) in the partial wave fit of ip{2S) 7XcO, XcO TT^Tr~ K~^K~ . 

The fitted mass and width of the /n(1710) are M = 1760 ± 15 MeV/c^ and F = 125 ± 25 MeV/c^. The mass is 
somewhat higher than the PDG value 0. The /o(1790) has been recently claimed in J/ip decays The parameters 
of /o(1790) and those of /o(1710) given by the PDG and Ref. [ill are tested in this analysis, and the log likelihood 
becomes worse by 10.3, 11.5, and 4.2, respectively. The log likehhood can be improved by 5.7 if both the /o(1790) and 
/o(1710) (using parameters of Ref. [ll|) are included into the fit, while the results for the other components only change 
a little. A fit replacing the decay modes with an /o(1710) with /2(1710) ~ namely, using cr/2(1710), /o(980)/2(1710), 
and /o(1370)/2(1710) decays instead of cr/o(1710), /o(980)/o(1710), and /o(1370)/o(1710) in the fit - is made in 
order to check the spin-parity of the structure around 1.75 GeV/c^ in the K^K~ mass region. It gives a worse log 
likelihood by 55.6 and a poor fit result, as shown in Fig. 7(a). 

There are two p-like resonances, /9(1450) and p(1700), which decay to KK final states in the 1600-MeV/c^ region 
As a check, we add p(770)p(1450) and p(770)/9(1700) intermediate decay modes into the fit with four more fitted 
parameters, using the mass and width measurements given by Ref. [H, and the log likelihood is improved by 6.4. 
Replacing the components cr/o(1710), /o(980)/o(1710), and /o(1370)/o(1710) with p(770)p(1450) and p(770)/9(1700) 
gives a worse log likelihood by 81.6, and the corresponding fit is shown in Fig. 7(b). An investigation of the p(770) 
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FIG. 7: Mass projections in K+ K' for fits replacing /o(1710) with (a) /2(1710) or (b) p(1450) and /9(1700). The dots with 
error bars are data, and the histograms are the fits. 



band in Fig. 4(a) indicates that the structures at high K^K mass are a reflection of the other intermediate decay 
modes. 

A free fit to /o(2200) gives a fitted mass of 2170 ± 20 McV/c^ and a width of 220 ± 60 MeV/c^. There are several 
scalar and tensor resonances located in the mass range greater than 2.0 GeV/c^ which decay to the K~^K~ final 
states Q: /3(2150), /2(2150), /o(2200), and /2(2300). A variety of alternative fits to the high K+K' mass region, 
listed in Table II, using these resonances and the /o(2100) have been tried. Note that PDG mass and width values 
are used in the fits for the /o(2100), /2(2150), /3(2150), and /2(2300) and that /o(2100), /2(2150), /o(2200), and 
/2(2300) couple with the a and /o(980), while the p(2150) couples with the p(770). We also tried /2(2200) instead of 
the /o(2200) in one fit, assuming it couples with the a and /o(980). 



TABLE II: Alternative fits to the high K+K' mass region instead of using only a /o(2200) involved in the fit; the right-hand 
column shows values of AS; positive values indicate poorer fits. 



Components 


AS 


a) 


p(2150) 


+35.3 


b) 


/2(2200) M=2170MeV/c^ 
r = 220 MeV/c^ 


+33.1 






c) 


/2 (2150) + /o (2200) 


-4.5 


d) 


/2(2150) -f /o(2100) 


+9.9 


e) 


/2 (2150) -f/2 (2300) 


+28.6 


f) 


/o(2100)-f /o(2200) 


-4.3 


g) 


/o(2100) -f/2 (2300) 


-3.5 


h) 


/o (2200) +/2 (2300) 


-9.0 


i) 


,f2(2150) 


+32.5 



In Table II, fit a), using p(770)p(2150) instead of cr/2(2200) and /o(980)/o(2200), gives a worse log likehhood by 
35.3. A fit with both /o(2200) and p{215Q) improves the log likelihood by about 3.0 compared to the case with 
/o(2200) only. Also b), e), and i), where no scalar is included, as well as d), give bad fits. Fits c), f), g), and h) 
improve the log likelihood a little, but not significantly. 

Based on these fit results, we conclude that a scalar state /o(2200) at about 2.2 GeV/c^ which decays to K'^K' 
is needed in this channel, but no additional significant 0"*""^ or 2"*""^ is required in this mass region according to our 
study in Table II. 



B. xco ^ {K+T^-){K--K+) 
The shape of the K*{%'d2) is described by a P-wave relativistic Breit-Wigner curve, with a width 



r^Fo 



mp 1 + r^pf^ 
m 1 + r^p^ 



P_ 

.Po_ 
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where m is the mass of the Kt: system, p is the momentum of kaon in the Kit system, Fq is the width of the resonance, 

1 Tp" 

niQ is the mass of the resonance, po is p evaluated at the resonance mass, r is the interaction radius and r— 2- 

1 + 

represents the contribution of the barrier factor. The value (3.4 ± 0.6 ± 0.3) (GeV/c)~^ measured by the K~tt~^ 
scattering experiment fl9l | as an approximate estimation of the interaction radius r is used. 

Adding a /C*(892)°i?*(l680)° decay mode improves log likchhood by 21.0, where the mass and width of K*{1680)° 
arc fixed to PDG values and the interaction radius r in the fit is set to the value 2.0 given by Ref. [l^ . 

Most of the peak at 1430 MeVyc^ is fitted as /sTq* (1430) A'j;; (1430) and Xq* (1430) A'2* (1430) + c.c, with only a 
smaU contribution from K^{U30)K^{U30) - the fitted number of A:|(1430)fs:2 (1430) events is about 11% of that of 
ATq (1430)^0 (1430). Adding the K^{U30)K^{U30) mode improves the log hkelihood by 7.1, which corresponds to 
about a significance of 3.3(t for two free parameters of the fitted amplitudes. The mass and width of the 1^2(1430) 
are fixed to the PDG values Q. A free fit to K^{U30) gives a fitted mass of 1455 ± 20 MeV/c^ and a width of 
270 ±45 MeV/c^. 

There is another structure visible around 1.95 GeV/c'^ in the Kir mass distribution. Among the three generalized 
C-parity allowed intermediate decays from Xco, i^o (1430)/?o(1950), K^{U30)K^{1950), and kK^{1950), the first 
decay gives a better log likelihood value in our fit than the other two by 22.9 and 14.6, respectively. A free fit gives 
a fitted mass of 1945 ± 30 MeV/c^, but the width is poorly determined, ~ 500 MeV/c^, so we adopt the PDG mass 
and width values for the _fsrQ(1950) in the partial wave analysis. 

In Fig. 4(b) there is a possible accumulation outside the (892)° A"* (892)° cluster. Its explanation may be the 
broad S-wave k. Since the properties of k are still controversial (cf., the review paper in the PDG Q, and references 
therein); we use a Breit-Wigner amplitude of constant width, without any phase space factor with the parameterization 
obtained in Ref. M= 790 MeV/c^, F = 860 McV/c^ to describe this broad structure. The log hkelihood becomes 
worse by 39.2 after removing the kk from the fit, and the corresponding projection in Kn mass, shown in Fig. 8, 
obviously disagrees with data at low Kir mass. 




FIG. 8: The Ktt mass projection for the fit without decay. 

Finally, an intermediate decay mode AT* (892)°/?* (2300)° is added into the fit. Here the AT* (2300)° is used as an 
effective structure to describe the poorly known high Ktt mass region. The log likelihood changes by 45.2 if the 
A:*(892)°i?*(2300)° is omitted. The scanned mass and width for A'* (2300)" arc - 2.3 GeV/c^ and 300 MeV/c^ 
respectively. No higher spin-parity tests such as 3~ or 5~ are made. 

C. XcO (A±7r+7r-)A=F 

In Fig. 5 there are strong A"i(1270) signals in both Kp{770) and Ar*(892)7r, and there is also a weak peak around 
1.4 GeV/c^ in the A:*(892)7r. In the partial wave analysis, A:i(1270)A: and A:i(1400)A: decays, where A:i(1270) 
decays to Kp{770), A:*(892)7r, and K^{U30)tt and A:i(1400) decays to A'*(892)7r, are added. The fit also requires an 
additional structure to describe the ATtttt mass distribution at about 1.4 GeV/c^. Here adding Ar(1460) in the fit will 
improve the log likelihood by 39.3; in the fit, the mass and width values given by the PDG [l| are used. The masses 
and widths of the A'i(1270) and A'i(1400) are also fixed to PDG values so as to reduce uncertainties. 

There are two lowest-lying axial-vector-meson octets. These correspond to the singlet (^Pi) and triplet (^Pi) spin 
configurations of two quarks in a P-wave orbital angular momentum state. The nonstrange, isospin 1 = 1 members 
of the two octets have opposite G parity: the &i(1235) is in the ^Pi octet and has G = +1, while the ai(1260) is in 
the ■^Pi octet and has G = —1. The strange members of the ^Pi and ^Pi octets, the Ka and Kb, respectively, are 
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mixtures of the observed physical states, the _R'i(1270) and the i^i(1400), where 

Ka = cos 6'A'i( 1400) + sin 6'/vi (1270), 
Kb = cos 61/^1 (1270) - sin gJ-iTi (1400), 

and the mixing angle is near 6* ~ 45° The dominant i4:i(1270) decay mode is to Kp{nQ) [B = 42% ± 6%); 

the A'i(1400) decays almost always to A'*(892)7r {B = 94% ± 6%). In the hmit of strict flavor SU(3) symmetry, the 
amplitudes for two-body decays to conjugate mesons in the same pair of octets should be equal. Thus, since decays 
to 6i7r are forbidden by G parity, decays to KbK are disallowed by SU(3), and one expects relatively pure KaK final 
states in Xco decays. And, since 9 ~ 45°, there should be roughly equal amounts of -ftri(1270) and i^i(1400). 

The remarkable feature of the Ktitt distribution is that the contribution of the ii'i(1400) in the final fit is very small 
compared to the large i4ri(1270) signal. Adding a Ki{lAQ{)) improves the log likelihood by 6.9 which corresponds to 
the significance of about 2.7cr. The partial wave analysis yields a ii'i(1270) KpiJlQ) signal of 68.3 ± 13.4 events 
and a i<:i(1400) A'*(892)7r signal of 19.7 ± 8.9 events. The function used to fit the 7^1(1270) is a convolution 
of an S-wave Breit-Wigner function for the i^i(1270) with a P-wave Breit-Wigner function for the /o(770) meson. 
Flavor-SU(3)-violating A'i(1270) — Ari(1400) asymmetry for XcO decay is observed. 

D. GOODNESS-OF-FIT 

To determine the goodness of fit, a is calculated by comparing histograms, e.g., a vector of Poisson distributed 
numbers n = (ni, ...^un), with a hypothesis for their expectation values Vi = E[ni]. As the distribution is Poisson 
with variances af ~ ui, the becomes Pearson' s statistic, 

2 _ (n-j - i^i) 



If the hypothesis v ~ (i^i, ...vn) is correct and if the measured values rii are sufficiently large, then the statistic 
will follow the x^ probability density function with the number of degrees of freedom (ndf) equal to the number of 
measurements N minus the number of fitted parameters. 

For an n-body final state, the number of independent kinematic variables is 3n — 4. Thus, one can compare 3n — 4 
selected independent variables with the fit results by defining a quantity 

3n-4 
2 _ \ ^ 2 

^all / J 

i=i 

to check the quality of a global fit, which obeys the x^ distribution approximately with the number of degrees of 
freedom equal to the total number of measurements minus the number of fitted parameters; and the individual x| 
will give a qualitative measure of the goodness of fit for each kinematic variable. 

The number of independent kinematic variables for ^^{23) 7XcOi XcO ~^ 'k~^it~ K~ process is 10 after the use of 
the 5-C fit with the additional XcO mass constraint. We use the following 10 distributions to check the goodness-of-fit: 

- {9, (t>)TT-KKK = the polar angle and azimuthal angle of the (tt+tt" AT+AT") in the laboratory system, 

- {9, (j))Tr-KK = the polar angle and azimuthal angle of the (7r+7r~Ar~) in the (tt+tt" A'+AT") center of mass, 

- (6*, (/))tx = the polar angle and azimuthal angle of the (tt+A'^) in the (tt+tt^A'^) center of mass, 

- [9, (/))^ = the polar angle and azimuthal angle of the 7r+ in the (tt+AT") center of mass, 

and the invariant masses for the (7r+7r~Ar~) and (tt+A'") systems. Table III shows the results, where the number of 
bins is taken as the number of degrees of freedom for each individual distribution, and Fig. 9 shows the projections 
of the 10 variables. There is excellent agreement between the data and fit. 

Using the results in Table III, the Xaj; obtained is 219.8 with 185 degrees of freedom (=237-52, where 237 is a sum 
of the bin numbers for the 10 distributions and 52 is the number of fitted amplitude parameters), which corresponds 
to a confidence level of 4%. 



VI. SYSTEMATIC ERROR 

In this analysis, the systematic errors are estimated by considering the following sources: 
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TABLE III: Check of goodness of fit using 10 independent kinematic variables, wliere ndf and C.L. are the number of degrees 
of freedom and the corresponding confidence level. 



Variable 


•J 

X 


ndf 


X /ndf 


C.L. 


COSO^^KK 


24.46 


18 


1.36 


0.14 


J. 


1 A OO 


on 


A ^1 
U. ( i 


U.OZ 


COsOttttK 


11.29 


20 


0.56 


0.94 


4>mrK 


22.58 


20 


1.13 


0.31 


cosOt^k 


11.34 


16 


0.71 


0.79 




18.81 


20 


0.94 


0.53 




12.15 


20 


0.61 


0.91 


4>TT 


12.63 


20 


0.63 


0.89 




37.26 


36 


1.04 


0.41 




55.07 


47 


1.17 


0.20 




(1) Uncertainty of the parameterization of the a line shape: use a Breit-Wigner amplitude of constant width M = 
470 MeV/c^, T = 613 MeV/c^ [13 instead of the Adler zero parameterization. 

(2) Uncertainty of the parameterization of the k line shape: c hang e the mass and width of the Breit-Wigner 
amplitude of constant width to M = 745 MeV/c^, T = 622 MeV/c^ [IJ 

(3) The parameters of the /o(980) are still uncertain, and in addition to the solution of Ref. '^^ ^^^o use 
measurements of some recent experiments such as E791, GAMS, and WA102 [2^ l23l |23 |. where a Breit-Wigner 
description with the width varying from 44 to 80 MeV/c^ was used for the /o(980). We determine the change both 
by using the solutions of Refs. [llllllll and by varying gi in Ref. [l3| from 0.1108 GeV/c^ to 0.090 GeV/c^ and 
0.130 GeV/c^ while keeping the ratio (72/51 fixed. 

(4) As mentioned above, we use the measurement of Ref. El, (3.4 ±0.6 ±0.3) (GeV/c)-i for r in the P-wave 
relativistic Breit-Wigner parameterization. We also use r varied by one sigma to 2.73 (GeV/c)~^ and 4.07 (GeV/c)~^ 
to determine the change in the fit. 

(5) Vary mass and width values of the (Kn) components /v^(1430), i4:*(1680), and /Vj;(1950) within the PDG 
errors 

(6) Vary mass and width values of the {Ktttt) components ^^1(1270), i^i(1400), and /\(1460) within the PDG 
errors Q, where the width of KiUGO) is changed to 200 and 300 MeV/c^ and mass to 1.40 and 1.46 GeV/c^ 

(7) Remove the small component /o(1710)/o(1370) from the fit, where /o(1710) decays to tt+tt" and /o(1370) decays 
to K+K-. 

(8) Add three {Ktt-k) resonances ^1(1650), ii:2(1770), and ii'2(1820) into the fit. 

(9) Remove the A'* (892)° A'* (2300)° + c.c. decay mode from the original fit but keep (8). Tests (8) and (9) strongly 
affect the branching fraction measurements. 

(10) Try alternative starting conditions for the maximum likelihood fit. 

(11) Uncertainty of the background in the partial wave analysis fit. In order to investigate the effect of the amount 
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and shape of background on the results, we increased the number of background events in the fit to 55, which is the 
maximum estimation, fitted without background, and added an incoherent phase space background normaUzed to 29 
events instead of fitting with the opposite sign of log likelihood for the MC background. The change on the fit results 
is very small, less than 5% for all the measurements. 

Total errors are obtained by adding the individual errors in quadrature. For the branching fraction uncertainties, 
the uncertainties in MDC tracking, kinematic fitting, PID, efhciency of the photon ID, and the number of tp{2S) 
events [2^ are also included, and the total systematic error for this part, 12%, is taken from Ref. 

VII. DISCUSSION 

From the XcO ^ tt^tt~ K~ decay fit results, it is found that scalar resonances have large decay fractions compared 
to those of tensors, and such decays provide a relatively clean laboratory in which to study the properties of the scalars 
/o(980), /p(1370), /o(1710), and /o(2200). There is also conspicuous production due to the X*(892)°A'*(892)" and 
A'o*(1430)i^^(1430) (and ii:*(1430)A'*(1430) + C.C.) pairs, and flavor-SU(3)-violating ii:i(1270) - A'i(1400) asymmetry 
is observed. Information on these states is very desirable and will be described below. 

For other components in Table I, because of the low statistics (/2(1270)/2(1270), /o(1710)/o(1370), 
and K 2 (14:30) K 2 (1^30)), uncertain parameters of intermediate resonances involved ((t/o(1710), ct/o(2200), 
A'* (892)" A'* (1680)° + c.c, ATq* (1430) A'^; (1950) + c.c, nk, and A'(1460)A'), and the poorly known high mass ATtt 
region which is described by the AT* (892)" AT* (2300)° + c.c, it is difficult to obtain precise quantitative results or 
make definite conclusions. The systematic errors on the numbers of events for decay modes /o(980)/o(1710) and 
/o(1370)/o(980) are very large, while the significance of the /o(980)/o(1710) is small; therefore we also will not 
consider these two decay modes in the following branching fraction measurements. 

The numbers used in the branching fraction (or upper limit) calculations and the corresponding results are 
summarized in Table IV, where the first errors are statistical and the second are systematic. The value of 
B['ip{2S) -> jXco], (9.22 ± 0.11 ± 0.46)%, recently measured by the CLEO experiment, is used |2||. The fit re- 
moving the Ari(1270) — > Kp{770) gives a worse log likelihood by 39.5, which corresponds to the signal significance of 
about 8.6(7. The 90% confidence level (C.L.) upper limit for XcO A'i(1400)A' is obtained by increasing the number 
of events by 1.28(t, where a includes the statistical and the systematic errors added in quadrature. 

TABLE IV: Summary of numbers used in the branching fraction (or upper limit) calculations and corresponding results, where 
X represents the intermediate decay modes, A'^-'^'' is the number of fitted events, and e is the efficiency. 



Decay mode 




ei%) 


Sys. error (%) 


B[xcO ^ X ^TT+n-K+K-] 
(xlO-*) 


Significance 


/o(980)/o(980) 


27.9 ±8.7 


6.25 ± 0.01 


+hh:/ 

-45.3 


3.46 ± l-OSti^^ 


5.3a 


/o(980)/o(2200) 


77.1 ± 13.0 


7.09 ± 0.01 


+ 19.6 
-27.2 


8.42 ± 1.42+2 29 


7.1a 


/o(1370)/o(1710) 


60.6 ±15.7 


6.59 ± 0.01 


+46.1 
-23.6 


7.12±1.85t?;g^ 


6.5a 


A:*(892)°A'*(892)° 


64.5 ± 13.5 


6.18 ±0.01 


+28.3 
-24.6 


8.09 ± 1.69+1-99 


7.1a 


Ao* (1430) A-Q* (1430) 


82.9 ± 18.8 


6.15 ±0.01 


+29.2 
-18.2 


10.44 ± 2.37t.i;gl 


7.2a 


Ao* (1430)^2* (1430) + c.c. 


62.0 ± 12.1 


5.66 ±0.01 


+15.6 
-23.4 


8.49 ± 1.66ti j9 


8.7a 


Ai (1270)+ A" + c.c. 












Ai(1270) ^ Ap(770) 


68.3 ± 13.4 


5.68 ±0.01 


+ 19.4 
-17.6 


9.32 ± 1.83li;^4 


8.6a 


Ai(1400)+A- +CC, 






+219 
-24.5 






Ai(1400) ^ A*(892)7r 


19.7 ±8.9 


4.94 ± 0.01 


< 11.9 (90% C.L.) 


2.7a 



The partial wave fit provides magnitudes and phases of the different partial amplitudes, as well as the interference 
terms. The intensity from these amplitudes is used to weight both the complete set of generated MC events and the 
set which survives the selection procedure. The ratio between these two weighted sets is the efficiency. 

Using the number of selected tp{2S) JXcO ^ Tjr+7r~ AT+AT" events, the overall efficiency determined by the 
method above, (5.85±0.01)%, and the result of Ref. [2g we get the corresponding branching fractions after subtracting 
background 

B[ip(2S) 7XcO ^ -/n+TT-K+R-] = (1.64 ± 0.05 ± 0.20) x 10'^, 
B[xco -K+Ti'K+K-] = (1.78 ± 0.05 ± 0.23) x 10"^ 

For the decay mode XcO ~* /o(980)/o(980), each /o(980) can decay to cither tt+tt^ or K^K~ , so it is necessary to 
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divide the result in Table IV by a factor of 2 to obtain the branching fractions 

B[i;{2S) ^ 7Xco ^ 7/o(980)/o(980)]S[/(980) ^ 7r+7r-]6[/(980) ^ K+ R-] = (1.59 ± 0.50t°;?^) x 10-^ 
S[xco /o(980)/o(980)]K[/(980) ^ ^+7r-]6[/(980) ^ - (1.73 ± QMt^^f^) x IQ-^. 

Combining this result and that of Ref. , we can determine the ratio of the partial decay width of /o (980) to tttt 
to those of TTTT and KK: 

r.. (6.5±1.93)x| 



^^■^ + ^Kk (6.5±1.93) X | + (1.59™) X 2 ' 

The numerical factors | and 2 take into account that (a) two-thirds of tttt decays are to tt+tt" and one-third to tt^tt^ 
and (b) there are equal numbers of decays to K^K~ and K'^K^ . Here the errors are the statistical and the systematic 
errors added in quadrature, and for the systematic error, the common parts related to the MDC tracking, kinematic 
fitting, efficiency of the photon ID, and the number of ip{2S) events cancel. 

There is a strong /o(980)/o(2200) with a signal significance of 7.1a in the XcO tt^tt~ K'^ K~ decay. The mass 
and width of the /o(2200) are optimized as M = 2170± 20ti^ MeV/c^ and T = 220 ± 60^4^ MeV/c^, and 

Bixco ^ /o(980)/o(2200)]6[/o(980) ^ 7r+^-]6[/o(2200) ^ K+R-] = (8.42 ± 1.42ti-|5) x 10"^ 

Changing the spin-parity of the /o(2200) in the fit or adding an additional resonance, shows that the spin-parity of the 
/o(2200) is well determined and no additional resonance is needed in the /o(2200) mass region. However, compared 
to the nearby states /o(2100) and /2(2150), its properties are still less well known P|, and more experimental data 
are needed. 

Another significant decay mode to /o(1370)/o(1710) is also foimd with a significance of 6.5(7, where /o(1370) decays 
to TT+TT" and /o(1710) decays to K+K'. The fitted mass and width of the /o(1710) are M = 1760 ± 15t}^ MeV/c^ 
and r = 125 ± 25^1° MeV/c^. The fitted mass is somewhat higher than the PDG value |]J. The spin component 
can be separated from spin 2 clearly, and replacing /o(1710) with either of the tensors /9(1450) or /9(1700) does not 
give a reasonable fit to the data. A free fit to /o(1370) gives a fitted mass of 1265 ± 30t^° MeV/c^ and a width of 
350 ± lOOlgQ^ MeV/c^. The corresponding branching fraction is 

B[xco -> /o(1370)/o(1710)]S[/o(1370) ^ ^+7r-]B[/o(1710) ^ K+K^] = (7.12 ± 1.85t?:^^) x 10^^- 

Besides the intermediate (7r+7r~)(iir+i^~) decay modes listed in Table I, we tried the following combinations in the 
fit: /o(1370)/o(1370), /o(1370)/o(1500), /o ( 1 500) /o( 1370), /o(1500)/o(1500), and /o(1500)/o(1710). None of them 
improved the log likelihood more than 5. So we didn't include these processes in the final solution of our fit and set 
upper limits at the 90% C.L.: 

S[Xco -> /o(1370)/o(1370)]S[/o(1370) ^ ^+7r-]S[/ol370) ^ if+if"] < 2.9 x 10-^ 
S[Xco ^ /o(1370)/o(1500)]S[/o(1370) ^ ^+7r]S[/o(1500) ^ K^K'] < 1.8 x 10-^ 
BiXcO -> /o(1500)/o(1370)]S[/o(1500) ^ ^+7r-]S[/ol370) ^ K+R-] < 1.4 x 10-^ 
B[XcQ ^ /o(1500)/o(1500)]6[/o(1500) ^ ^+7r-]6[/o(1500) ^ K+R-] < 0.55 x IQ-^, 
B[xco ^ /o(1500)/o(1710)]S[/o(1500) ^ ^+7r-]6[/o(1710) ^ K+R-] < 0.73 x 10"^ 

From the results of the fit, the branching fraction B[xcO ^ ii:*(892)° A'*(892)° tt+tt~K+K-] ^ (8.09 ± 
1.69lf 

'99) X 10 ^ is obtained. Using the branching fraction of K*{892)'^ to the charged Kn mode, which is taken as 
I, we get 



B[^|:{2S) jXcO -fK*{892fK*{892)"] = (1.68 ± 0.35t^;4^) x 10 

1.52 
1.45 



B[xco K*{892fK*{892f] = (1.82 ± O.SSto.ll) x 10"^ 



The values obtained here are consistent with those in Ref. 

Most of the peak around a Ktt mass - 1430 MeV/c^ is fitted with iro*(1430)^^(1430) and i4ro*(1430)if|(1430) 
-f- c.c, with only a small contribution from A"|( 1430) A'2 (1430). The measured branching fractions are 

B[XcO A:o(1430)A'fJ(1430)° ^ tt+tt- A'+A^] = (10.44 ± 2.37+i^;°g) x 10-^ 
B[xcO -> A"* (1430) A'2 (1430)° + c.c. ^ tt+tt-K+R-] = (8.49 ± IMtlH) x 10""*. 
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A free fit to A"* (1430) gives a fitted mass of 1455 ± 20 ± 15 MeV/c^ and a width of 270 ± 45^^^ MeV/c^. The 
fitted number of /q (1430)^^2(1430) events is about 11% of that of i^g (1430)^o (1430), and the signal significance of 
i4r|(1430)i^2(1430) mode is about 3.3tT. These measurements are important for the study of Xco decays, as well as 
those of K*(1430) and i^2*(1430). 

Using the results shown in Table IV and the branching fractions of A'i(1270) Kp{770) and ^1(1400) K*{892)n 
given by the PDG [2, we determine 

BiXco Ki{l270)+K- + c.c] = (6.66 ± 1.3ll} j?) x 10-^ 
B[xcQ Ki{UOO)+K- + c.c] < 2.85 x 10-^ 

at the 90% C.L.. To accommodate this, a mixing angle of > 57° is required. A flavor-SU(3)-violating 7^1(1270) — 
i4ri(1400) asymmetry is observed. The asymmetries with opposite character for ■0(25') and J/ip decays were also 
observed in the BESI data where for ^(25') data, 9 < 29° and for J/tp data 9 > 48°. 

VIII. SUMMARY 

In summary, a partial wave analysis on XcO ~^ tt^tt^ K~ in ^p{2S) — > 7XcO decay is performed using a sample of 14 
million '0(25') events. From the fit we find significant contributions to the XcO decays from the channels /o(980)/o(980), 
/o(980)/o(2200), /o(1370)/o(1710), K*{892fK*{892f, K^{U30)K^{U30), i^^(1430)i4:2*(1430)+c.c. and Ki{1270)K. 
The mass and width of the /o(1710) are determined to be 1760 ± 15t\l MeV/c^ and 125 ± 25^11 MeV/c^, and those 
of the /o(2200) are 2170 ± 20^1^ MeV/c^ and 220 ± 60^4^ MeV/c^. Flavor-SU(3)-violating A'i(1270) - A'i(1400) 
asymmetry is observed, with the mixing angle 9 > 57°. 
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